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A solution containing ^^Ge in equilibrium 
with its daughter, ^^Ga, has been 
standardized for the first time at the 
National Institute of Standards and 
Technology (NIST) using 3 hquid 
scintillation-based techniques: live-timed 
A%p -y anticoincidence (LTAC) counting, 
the Triple-to-Double Coincidence Ratio 
(TDCR) method, and ^H-standard 
efficiency tracing with the CEEMAT / 
NIST (CNET) method. The LTAC 
technique is much less dependent on level 
scheme data and model-dependent 
parameters and was thus able to provide a 
reference activity concentration value for 
the master solution with a combined 
standard uncertainty of about 0.3 %. The 
other two methods gave activity 
concentration values with respective 
differences from the reference value 
of +1.2 % and -1.5 %, which were still 
within the experimental uncertainties. 



Measurements made on the NIST "47i"7 
secondary standard ionization chamber 
allowed for the determination of calibration 
factors for that instrument, allowing future 
calibrations to be made for ^^Ge/^^Ga with- 
out the need for a primary measurement. 
The ability to produce standardized 
solutions of ^^Ge presents opportunities 
for the development of a number of NIST- 
traceable calibration sources with very low 
(<1 %) relative standard uncertainties that 
can be used in diagnostic medical imaging. 
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1. Introduction 

The use of Positron Emission Tomography (PET) as 
a tool for diagnosing diseases, particularly cancer, 
continues to rise at a rate of about 20 % per year [1], 
partially due to technological advances that allow for 
more quantitative data to be collected. The ability to 
consistently acquire truly quantitative imaging data 
depends on the use of radioactivity standards traceable 
to the National Institute of Standards and Technology 
(NIST). 

Currently the most widely used radionuclide in 
PET imaging is ^^F. Although NIST has previously 



^ CIEMAT is an acronym for Centro de Investigaciones Energeticas, 
Medioambientales y Tecnologicas, which is the National Metrology 
Institute of Spain. 



standardized ^^F [2,3], its short half-hfe (1.83 h) makes 
it nearly impossible to prepare and distribute Standard 
Reference Materials to most users. As a potential solu- 
tion, ^^Ge, in equilibrium with its daughter, ^^Ga, has 
been proposed as a surrogate. 

A simplified scheme for the decays of ^^Ge 
and its ^^Ga daughter is shown in Fig. 1. Germanium-68 
decays by pure electron capture (EC) to the ground 
state of ^^Ga with a half-life of 270.95(16) d 
[4]. Gallium-68 in turn decays with a half-life of 
67.71(9) min by a combination of EC and positron 
emission primarily to the ground state of ^^Zn, but also 
with a branch to an excited state at 1077 keV with 
a probability of about 3 % and a number of higher 
excited states with a combined probability of 
under 0.4 %. 
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Fig. 1. Simplified decay scheme for Ge- Ga. Data were taken from the DDEP data evaluation [4]. The 
branching probabilities b^, b2, and b3 refer to the two positron emission and gamma emission probabili- 
ties, respectively. 



The nature of the decay scheme of ^^Ge/^^Ga makes 
it amenable to a variety of different standardization 
techniques. In 1994, Schonfeld, et al. [5] reported on 
the results of measurements made with liquid scintil- 
lation (LS) counting using the CIEMAT/NIST ^H 
standard efficiency tracing method (CNET) [6,7], 
Anfi - 7 coincidence, and a calibrated ionization 
chamber (IC). The data showed good agreement 
between all three methods, giving activity concentra- 
tion values within the respective experimental uncer- 
tainties (nominally 1 % relative standard uncertainty). 

More recently, Grigorescu, et al. [8], reported on the 
results of measurements using Anfi - y coincidence 
counting. As with Schonfeld, et al., the coincidence 
spectrometer consisted of a proportional counter and 
Nal(Tl) detector for the ^ and y detection channels, 
respectively. Because this experimental arrangement 
requires the use of dried sources, corrections for the 
loss of ^^Ge due to chemical volatility were necessary in 
both studies. This effect is reported by Grigorescu to be 
on the order of 20 % to 26 %. Nonetheless, they were 
able to obtain a measurement result with about a 1 % 
relative standard uncertainty. 

Liquid scintillation (LS) counting has been the 
method of choice in our laboratory for the measurement 
of j8-emitting radionuclides, primarily due to the high 
LS detection efficiency and the relative ease of sample 
preparation. Methods based on LS counting have 
another advantage in the context of measuring ^^Ge 
because the sample is introduced into the cocktail while 



still in solution, thereby eliminating the need to prepare 
dry sources. Seeking to take advantage of this, we have 
measured a single solution of ^^Ge/^^Ga using three LS- 
based methods: live-timed 471^8 - 7 anticoincidence 
(LTAC) using LS as the j8 counting channel, the Triple- 
to-Double Coincidence Ratio Method (TDCR) [9,10], 
and CNET [6,7]. 



2, Experimental 

All evaluation of measurement uncertainties 
throughout this work follow accepted conventions used 
by the NIST Radioactivity Group and are in accordance 
with those recommended by the principal metrology 
organizations [11]. All individual uncertainty compo- 
nents are given as estimated experimental standard 
deviations (or standard deviations of the mean, if 
appropriate), or quantities assumed to correspond to 
standard deviations regardless of the method used to 
evaluate their magnitude. Unless explicitly stated, all 
uncertainties cited in this paper are "standard uncertain- 
ties," corresponding to one uncertainty interval. One 
particular exception is the uncertainty reported for the 
activity concentration of the calibrated ^^Ge solution, 
which is given as an "expanded combined standard 
uncertainty." In accordance with NIST pohcy [12], the 
combined standard uncertainty (calculated by combin- 
ing the individual uncertainty components in quadra- 
ture) is multiplied by a "coverage factor" of ^ = 2 to 



266 



Volume 1 13, Number 5, September-October 2008 

Journal of Research of the National Institute of Standards and Technology 



obtain an "expanded uncertainty" assumed to give an 
uncertainty interval having a confidence level of 90 % 
to 95 %. 

2.1 Initial Solution Preparation 

The master solution used in these experiments 
contained nominally 125 MBq ^^Ge in 5 mL of 
0.5 mol ■ L"^ HCl and was prepared by International 
Isotopes Idaho, Inc. (Idaho Falls, ID)^ using ^^Ge 
produced at the 1 00 MeV Isotope Production Facility 
at Los Alamos National Laboratory using the 
""'Ga(p,2n)^^Ge reaction. 

A generalized scheme for the preparation of the 
counting samples is shown in Fig. 2. The first step 
involved transfer of the master solution out of the ship- 
ping vial into a NIST standard 5 mL flame-sealed 
ampoule while at the same time performing the first of 
three serial dilutions that would be needed in order to 
bring the activity level in one of the ampoules down to 
that suitable for LS counting. The ampoule used for this 
study labeled A], was prepared by volumetrically 
adding 1 mL of the stock solution to an ampoule con- 
taining 4 g of gravimetrically added carrier solution 
having nominally 45 |ng each of nonradioactive Ge^ 
and Ga^^ per gram of solution using 0.5 mol • L"^ HCl 
as the solvent. The ampoule was weighed again after 
the addition of the ^^Ge to determine the mass of added 
radioactive solution. Ampoule A2, shown in Fig. 2, was 
held in reserve for future experiments. 

Ampoule AJ was measured in the NIST-maintained 
radionuclide activity calibrators ("dose calibrators") 
and the NIST "47t"7 IC [13] to allow for the determina- 
tion of calibration factors in this specific measurement 
geometry. The solution in ^7 was then diluted by a 
factor of about 200 through two serial gravimetric dilu- 
tions, giving two additional ampoules, AJDJ and 
A1D2. As an additional check of the dilution factor 
between A J and AJDJ, the latter was also measured in 
the NIST IC. 

2.2 Liquid Scintillation Source Preparation 

All counting sources for the three counting tech- 
niques were prepared using solution AJD2. A total of 
1 8 LS cocktails containing ^^Ge were prepared for these 
studies. For the LTAC experiments, 4 mL of HiSafe-3 
(Perkin Elmer, Waltham, MA) or PCS (GE Healthcare 



Certain commercial equipment, instruments, or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the purpose. 



Biosciences, Piscataway, NJ) scintillant were added to 
each of two 3 cm diameter, glass pseudo-hemispheres. 
Nominally 0.04 g of ^^Ge solution were gravimetrically 
added and the hemispheres were sealed using epoxy 

For the TDCR experiments, two cocktails each of 
HiSafe-3 and PCS were prepared by dispensing 10 mL 
of the scintillant into four 22 mL borosilicate glass LS 
vials (two per scintillant), followed by the gravimetric 
addition of nominally 0.04 g of solution from AJD2. 
Similarly, 10 cocktails were prepared for the CNET 
experiments (five vials per scintillant). To vary the 
counting efficiency of the CNET cocktails, between 

2 drops and 18 drops of a 10:1 (by volume) dilution of 
nitromethane in ethanol were added as a quenching 
agent to the CNET cocktails in addition to the scintil- 
lant and radioactive solution. In order to perform the 
efficiency tracing, a separate set of 1 LS vials having 
identical composition to the ^^Ge cocktails were pre- 
pared using a dilution of a NIST tritiated water 
Standard Reference Material 4927F [14] in place of the 
^^Ge. In order to make the ^^Ge and "^H cocktails as close 
in composition as possible, nominally 1 mL of the 
Ge'^VGa'^^ carrier was added to each of the cocktails. 

Two background blanks (one for each scintillant) 
were prepared for the TDCR and LTAC measurements 
in their appropriate vials. In order to properly mimic 
the composition of the radioactive cocktails, an equiv- 
alent mass of nonradioactive carrier solution was added 
to each blank. For the CNET measurements, four 
blanks were prepared so as to have the identical sample 
compositions of the least- and most-quenched of the 
radioactive ^^Ge cocktails. As with the TDCR and 
LTAC blanks, nonradioactive Ge/Ga carrier was substi- 
tuted for the ^^Ge solution. 

2.3 4nP- /Anticoincidence Counting (LTAC) 

The system constructed at NIST uses an LS source 
optically coupled to an appropriate photomultiplier 
tube for the beta channel and a thallium-doped sodium 
iodide [Nal(Tl)] detector for the y-ray channel, as 
described previously [15,16]. The LS-based system is 
well suited to this case since the source solution does 
not have to be dried, and therefore, the large (20 % to 
26 %) correction for Ge loss reported by Grigorescu [8] 
is avoided. 

The four active sources were each measured for 
between 2 and 4 cycles and the blank sources up to 

3 cycles during the period from 24 April to 5 May 2007. 
Each counting cycle consisted of measurements at 
between 8 to 12 threshold levels on the LS detector 
for between 200 and 1000 seconds. A minimum of 
5-10^ LS and 1-10^ anti-coincident Nal detector 
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Fig. 2. Scheme used for preparation of sources used to calibrate solutions of ^^GeCl4. 
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counts were recorded for each non-blank measure- 
ment. The LS signal-to-background ratio for the lowest 
threshold (highest background) data points was about 
950:1, while the signal-to-background ratio for the 
Nal detector was about 350:1. Further systematic tests 
demonstrated that the background variability during 
the span of the experiment, the variation of extend- 
ing dead-time, and the presence or absence of the 
aluminum absorber did not affect the measurement 
results. 

The positron decay of the ^^Ga was detected in the 
LS channel, with count rate A^^ , while electron capture 
events from both ^^Ga and ^^Ge were avoided by con- 
straining the lower level discriminator (LLD) on the 
amplified signal to be above about 20 keV beta energy. 
In this way, the determined activity value was inde- 
pendent of any atomic transitions (all below 1 1 keV), 
and directly proportional to the total positron emission 
probability. The LS positron efficiency, £p, was varied 
between about 0.5 and 0.95 using the LLD and extrap- 
olated to 1.0. The Nal(Tl) detector was gated on the 
511 keV region using a single channel analyzer and the 
total 7-ray (A^^ ) and anticoincidence (N^c) count rates 
were recorded. The extrapolation parameter used was 
Y= N^c I ^y ^ (1 - ^^)- Most of the T^ray counts were 
due to positron-annihilation decays, detected with effi- 
ciency. 



N,-K{b,+b,)\\ + 



^71077 b^ ~ t>2 

^.n. b, +h 

aim I Z 



(3) 



A^. 



N,{b,+b,) 



0) 



where A^q is the activity and (Z?i + b^ is the total positron 
emission probability. There was an additional, approxi- 
mately 0.2 %, contribution from Compton scat- 
tering of 1077 keV 7-rays, detected with efficiency 
£^1077- Since some of these 1077 keV 7-rays correspond 
to electron capture events, and not positron emission, a 
small (0.2 %) correction to the intercept was necessary. 
The modified extrapolation formula is. 



Equation (3) is similar to Eq. (4) in [8], with a signif- 
icant difference that here we account for the fact that 
some of the 1 077 keV 7-rays do correspond to a 
positron branch, and thus do not need to be subtracted. 
This effect is accommodated by the presence of Z?2 in 
the numerator of the correction term (final term in 
Eq. (3)). Corrections due to the LS efficiency for 7-rays 
and for coincidences due to those events were not nec- 
essary, as described below. 

An estimate of £^1077 was obtained during a separate 
set of measurements with a ^^Co point-source and ^^m 
was estimated using Eq. (1). The correction term was 
checked by exaggerating the effect during additional 
measurements made with various 7-ray energy gates. 
Gates G7, G2, and G3 corresponded to a gate around 
the Compton region of the annihilation spectrum, the 
annihilation photopeak (511 keV) and the annihilation 
sum peak (1022 keV) respectively. The sum peak was 
unresolved from the 1077 keV peak in the Nal(Tl) 
detector events. As designed, the data from gates Gl 
and G3 required large corrections due to reduced ^^m 
and enhanced £^1077 respectively. The uncorrected and 
corrected A^q values are shown in Table 1 . The fact that 
the large corrections for Gl and G3 gave consistent 
results with G2, supports the use of this method for the 
small (0.2 %) correction to the final A^q value, based on 
G2 alone. 



Table 1, Calculated ^^Ga activity from various 7-ray gates and using 
Eq. (3), with and without the final correction term, and relative to the 
corrected G2 value, (used for the final activity determination). The 
uncertainties {k= 1) on the uncorrected activities are standard devia- 
tions of the intercepts from the least-squares fits to the various data 
sets. The uncertainties (^=1) of the corrected values are from 
estimates of the uncertainties in e^^^ and £^1077 
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y-ray gate 



Gl - Comptons 
G2 - photopeak 
G3 - sum peak 



Uncorrected A^o (Bq) Corrected A^o (Bq) 



1.020 ±0.002 

1.0021 ±0.001 
1.033 ±0.002 



1.002 ±0.004 
1.0000 ±0.0004 
1.000 ±0.007 



where the numerator and denominator of the correction 
term correspond to total 7-rays and anticoincident 
7-rays, respectively and the branching probabilities 
by, Z?2, and Z?3 are illustrated and enumerated in Fig. 1. 
Note that the extrapolation is linear in 7, and the 
7= (£^ = 1) intercept is given by 



No correction in the final result was made for the 
gamma efficiency of the LS detector, or for coinci- 
dences between such events and the Nal (Tl) detector. 
For gates Gl and G2, these effects tend to cancel out 
due to the two-photon annihilation process. If one 
photon is detected in the LS detector, that efficiency 
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can be monitored by the other photon interacting in the 
Nal detector [15]. If such an effect were present it 
would lead to a non-linear efficiency extrapolation, A 
typical G2 data set and residuals from a linear least- 
squares fit are shown in Fig, 3 and it is evident that a 
linear fit is satisfactory. Yet, a quadratic extrapolation 
was needed to fit the entire range of G3 (sum peak) data 
due to the unmonitored LS efficiency for 7-rays in 
that configuration. Thus, for G3, a smaller £^ range 
(0.9 - 0,95) was employed such that linear and quadrat- 
ic fits gave consistent results. This value was only used 
for the £yio77 sensitivity test. 

Another possible cause for a non-linear extrapolation 
would be if both the LS and Nal(Tl) efficiencies dif- 
fered depending on whether the positron was stopped in 
the LS hemisphere, or escaped before annihilating. This 
effect was mitigated by three factors: (1) the fact that 
most positrons annihilated within the hemisphere, (2) 
the high LS efficiency, and (3) the well-type geometry 
of the Nal(Tl) detector. The sensitivity of the result to 
this effect was tested by placing an approximately 
0.5 cm thick aluminum foil over the hemisphere 



and comparing the resulting activity determination. 
No change in the goodness of the linear fit was detect- 
ed and the ratio of the intercept with to without the 
foil was 1.000 ±0,001, where the uncertainty is 
a standard (^=1) uncertainty on the linear fit coeffi- 
cients. 

2.4 Liquid Scintillation Counting Using the Triple- 
to-Double Coincidence Ratio (TDCR) Method 

Each counting source was counted in the NIST 
TDCR system [17] on at least two separate occasions 
over the course of 27 days. Counting times were typi- 
cally 1200 s, which allowed for the accumulation of at 
least 1 0^ counts in each of the three doubles counting 
channels. For each counting experiment, data were 
acquired at a minimum of 4 efficiency points, which 
were varied through the use of a set of grey filters that 
were fitted over the LS vials. Data were acquired in 
triplicate at each efficiency point. The experimental 
efficiencies for the logical sum of double photon 
coincidence events, ^lsd? ranged from 0.89 to 1.14, 




0.6 0.7 



Fig. 3. Typical ^^Ga LS rate, A^^ , extrapolation versus anticoincidence efficiency, Y (bottom) and residual, /?, 
from a linear least-squares fit versus 7 (top). Statistical ik= i) uncertainties are shown in the residual plot. 
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The counting data were analyzed using a program 
developed in-house for use with the Mathematica [18] 
symbolic mathematics package. Details of the program 
and the computation strategy will be published sepa- 
rately However, it should be noted that the program 
calculates the total detection efficiency for the case of 
decay of ^^Ge in equilibrium with its ^^Ga daughter. To 
do this, the program was required to solve the TDCR 
equations [9,10,19,20] for the EC branch of the ^^Ge 
parent, as well as both the EC and positron decay of the 
^^Ga daughter. A relatively simplistic model, consider- 
ing twelve possible decay pathways, was adopted to 
describe the atomic transitions encountered in the EC 
decay of ^^Ge and ^^Ga. These are depicted in Fig. 4. 
The values of the various nuclear and atomic input data 



were taken from the evaluation of the Decay Data 
Evaluation Project (DDEP) [4]. 

The analysis program calculates the individual 
phototube efficiencies, thereby allowing for correction 
due to asymmetry in the counting rates in each of the- 
doubles counting channels. The contribution to the 
detection efficiency due to detection of the 511 keV 
anihillation photons was taken into account by using 
the positron spectrum calculated by the program 
SPEBETA [21] as input for the Monte Carlo simulation 
package PENELOPE [22] using the techniques 
described in [23]. The resulting spectrum of energy 
(positrons+annihilation photons) absorbed in the LS 
cocktail was then used as input data for the TDCR 
analysis code. 
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Fig. 4. Energy decay schemes used in calculating TDCR efficiencies for electron capture branches in the decay of ^^Ge and ^^Ga. 
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The stopping power, dE/dx, for electrons in the LS 
cocktail was calculated by fitting a function of the form 



(dE/dx) = a + bE + c(\nEf + d\nE/E + e/E (4) 



(E is the value of the midpoint energy for each bin of 
the calculated beta spectrum and a, b, c, d, and e are 
fitting parameters) to data from the NIST ESTAR [24] 
database using previously published LS cocktail 
compositions [25]. 

A separate program, assuming equal phototube 
efficiencies, was developed for evaluating the effects 
of varying different input and model parameters. 
Calculations of £lsd were made as a function of 
the TDCR for kB values between 0.009 cm ■ MeV"^ 
and 0.018 cm • Me V"^ and the resulting £lsd values 
were found to be insensitive to the value of kB. For 
consistency with previous measurements made in 
this laboratory [26], the value of kB for all analyses 
was taken to be 0.012 cm • MeV"\ A plot of the 
theoretical £ lsd values as a function of TDCR at 
kB = 0.012 cm ■ MeV"^ is shown in Fig. 5. 

2.5 Liquid Scintillation Counting Using the 
CIEMAT-NIST ^H-Standard Efficiency 
Tracing (CNET) Method 

Each LS cocktail was sequentially counted for 
10 cycles of 25 min per source in a Packard 
(Perkin Elmer, Waltham, MA) 2500TRLS spectro- 
meter. Samples were then removed from the counter, 
agitated and sequentially counted for 10 cycles of 
30 min per source in a Beckman LS6500 (Beckman 
Coulter, Fullerton, CA) spectrometer. 

Efficiency tracing involves calculating a relation- 
ship between the measured ^H LS efficiencies and 
the LS efficiencies expected for ^^Ge, in equilibrium 
with its daughter ^^Ga, over a range of experimental 
quench indicating parameters [7]. The efficiency 
tracing computer program CN2004 [27] was used in 
the analysis of the LS data after changing the default 
input file to include the nuclear and atomic data found 
in the DDEP evaluation [4]. The average calculated 
^^Ge/^^Ga efficiency was nominally 138% in the 
Packard LS counter and 147% in the Wallac LS 



counter using a kB value of 0.0075 cm • Mev"^ and 
assuming that the cocktail had the composition of Ultima 
Gold as specified in the default CN2004 input files. Plots 
of the calculated theoretical ^^Ge, ^^Ga, and total efficien- 
cies as a function of "^H tritium are given in Fig. 6. 

2.6 Ionization Chamber Measurements 

For the NIST IC measurements, both A J and 
AID] were measured 40 times each, in four groups 
of 10 measurements, alternating with 5 groups of 
10 measurements of either radium (^^^Ra) reference 
source RRSIOO or RRS500b. Results are analyzed as a 
ratio of the response of the ampoule to the response of 
the RRS. After correction for background, the resulting 
ratio is used to derive a calibration factor, or K-value, 
defined as the activity of a given radionuclide that 
would produce the same response as the RRS. The 
relative values of the RRSIOO and RRS500b are well 
characterized. By determining the K-value using the 
activity derived from different ampoules of different 
activity levels, it is also possible to verify the gravimet- 
ric dilution factor. The dilution factor from ^7 to AlDl 
was verified by this method to within 0.022 %>. The 
LTAC activity values and the mass dispensed into 
Al were used to determine K- values that can be used 
for future measurements of ^^Ge in the NIST ampoule 
geometry. 

2.7 Gamma Ray Spectrometry 

The solution that remained in AlDl after making the 
LS cocktails was analyzed for possible photon-emitting 
radionuclidic impurities using two calibrated High- 
Purity Germanium (HPGe) photon spectrometers at 
two different counting distances each. In addition, the 
data provided an additional, confirmatory measurement 
of the activity concentration using the 1078 keV 
gamma ray from the decay of ^^Ga. Characteristics of 
the detectors used in this study are given in Table 2. 

Data were collected using the GammaVision-32 
(Ortec, Oak Ridge, TN) software package and analyzed 
using both Gamma Vision-32 and Genie 2000 
(Canberra, Meriden, CT). Detection efficiencies were 
calculated from efficiency-energy relationships deter- 
mined using solutions previously calibrated at NIST 
and measured in the 5 mL NIST ampoule geometry. 
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Table 2. Characteristics of HPGe detectors used in the present study 



Detector parameter 



X-detector 



B-detector 



Detector diameter 

Detector length 

End cap window material 

Window thickness 

Crystal-window distance 

Crystal top dead zone thickness 

Crystal material 

Crystal hole depth 

Crystal hole diameter 

Detector side cap thickness 

Detector side cap diameter 

Detector side cap material 

Detector type 

Calibration Geometries (distances are source-to-detector) 



43.6 ±0.1 mm 

36.2 ±0.1 mm 

Beryllium 

0.5 ±0.05 mm 

3 ± 0.5 mm 

0.3 ± 0.03 [am 

Germanium 

32.6 ± 0.5 mm 

10.4 ±0.5 mm 

1.3 ±0.1 mm 

70 ± 1 mm 

Aluminum 

n-type 

Ampoule; side-mount, end-on 24 cm 



54.9 ±0.1 mm 

54.2 ±0.05 mm 

Beryllium 

0.5 ± 0.05 mm 

3 ±0.5 mm 

0.3 ± 0.03 [am 

Germanium 

47.2 ± 0.5 mm 

12 ±0.5 mm 

1.3 ±0.1 mm 

63.5 ± 0.5 mm 

Magnesium 

n-type 

Ampoule; 24 cm, 35 cm 



3. Results and Discussion 

3.1 Impurity Analyses 

No photon-emitting radionuclidic impurities were 
detected in solution ^7 to within the following limits 
(at the reference time) of the massic photon emission 
rate: 

785s-' -g-' for 30keV < E <507keV; and 
285s-' -g-' for 515keV<^<1800keV; 

where E is the gamma-ray energy. 

3.2 Activity Measurements Results 

The results of the massic activity determinations for 
the solution contained in ^7 as of the reference time 
are given in Table 3. The values in Table 3 take into 
account the dilution factor of 204.903231 between the 
solution in AI and that used in the assays, A1D2. The 
uncertainties given in the table are expanded (k = 2) 
uncertainties based on the components given in Tables 
4-7. 

Of the different techniques used in this study 
to determine the activity concentration of the ^^Ge 
solution, the LTAC technique is much less dependent 
on level scheme data and model parameters not direct- 
ly measured in the experiment. For this particular 
measurement, the only input parameter significantly 
impacting the activity calculation that was not directly 
measured in the experiment was the positron branching 
ratio. The other branching ratios only contributed to the 



Table 3. Results of massic activity determinations (C^, in Bq ■ g"^) 
for the ^^Ge solution contained in ampoule Al as of the reference 
time of 12:00 EST 1 May 2007. The uncertainties, given in parenthe- 
ses, are expanded (k = 2) uncertainties based on the evaluated uncer- 
tainty components listed in Tables 4-7 for the respective techniques 



Technique 



|6 T3^ . „-i 



Ca, lO^Bq-g 



4Kp - y anticoincidence counting 
(LTAC) 

LS counting with the Triple -to -Double 
Ratio (TDCR) method 



LS counting with the CIEMAT/NIST 
■^H-standard efficiencytracing method 
(CNET) 



Gamma-ray spectrophotometry with 
High Purity Germanium 
(HPGe) detectors 



3.104(18) 



3.141(25) 



3.058(44) 



3.2(9) 



minor (0.2 %) correction for the leakage of 1077 keV 
7-rays into the annihilation y-ray gate. And even this 
small correction was checked experimentally by modi- 
fying the experimental design to exaggerate the effect 
and then verifying that the corrected activity agreed 
with the original value. 

On the other hand, our implementations of the TDCR 
and CNET efficiency tracing methods are unable to 
separate the positron and EC decay signals and must 
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Table 4, Uncertainty components evaluated in the determination of the massic activity, C^, for ^^Ge solution ^7 by 4Kp - y anticoincidence 
counting (LTAC) 



Component, w^- 



Comment 



Evaluation type 



% 



Measurement variability 



Standard deviation of the mean on the 
determination of C^ for 8 trials 
encompassing 4 samples, 3 backgrounds, 
2 extending dead-times and 
absorber/no absorber- over a 
2-week period 



0.03 



Background variability 



Additional variation estimate based 
on 3 trials and within-run variation 



0.03 



Additional variability in 
background 



Additional variability not embodied 
by "random sources" 



0.05 



Half-life 



Standard uncertainty in half-life . 
(0.059 %) over the measurement 
decay interval 



0.0008 



Livetime 



Estimated standard uncertainty on C^ 
due to uncertainty in counting livetime. 



0.1 



Branching ratio 



Estimated standard uncertainty due to 
uncertainty in published decay 
branching ratios 



0.13 



Extrapolation 



Estimated standard uncertainty due to 
extrapolation to zero non-detection 
efficiency; based on sensitivity tests, 
previous measurements and models 



0.2 



Correction due to 
detection of 1077 keV 
photons 

Mass determinations 



Estimated standard uncertainty on C^ due 
to efficiency of detecting 1077 keV 
photons based on sensitivity tests 

Estimated standard uncertainty of mass 
for any single LS cocktail 



0.1 



0.05 



Dilution factor 



Uncertainty in C^ of solution GelAl due to 
uncertainty in gravimetrically-determined 
dilution factor between solutions in 
ampoules GelAl and GelAlD2 



0.04 



Combined {Uq = ^Ym]) 



0.29 



Expanded {Uq = Uq' k;k = 2) 



0.58 
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Table 5. Uncertainty components evaluated in the determination of the massic activity, C^, for ^^Ge solution^/ by hquid scintillation counting 
using TDCR method 



Component, Ui 



Comment 



Evaluation type 



% 



Sample repeatability 



Standard deviation of the mean on the 
determination of massic activity for a 
single LS cocktail (?7 = 3 - 5 
determinations of C^ per source) 



0.06 



LS cocktail composition 
variability 



Standard deviation on the determination 
of C^ for three LS cocktail compositions 
{n = 6-\9 determinations of C^ 
per composition) 



0.09 



Efficiency dependence 



Effect of ^^Ga beta 
endpoint energy, E^ ^^^ on 
eflSciency calculations 

Effect of other atomic and 
nuclear input data 



Half-hfe 



Mass determinations 



Median difference between maximum and 
minimum value of C^ determined for a single 
source at between 3 and 4 eflSciency values, 
varied by use of grey filters {n = 12 
independent measurements) 

Standard uncertainty in efficiency 
calculation due to standard uncertainties 
on positron endpoint energies of ^^Ga 

Standard uncertainty due to uncertainties on data 
used as input to the TDCR analysis code as 
determined by Monte Carlo methods. A total of 20 
data sets were generated from normal distributions 
defined by the pubhshed nuclear and atomic data 
and their associated standard uncertainties, which 
were taken as the standard deviation of the respective 
distributions. Each data set was used to calculate C^ 
using a single experimental data set 

Standard uncertainty in half-hfe (0.059 %) 
over the measurement decay interval 

Estunated standard uncertainty of mass 
for any single LS cocktail 



0.20 



0.18 



0.26 



6-10-3 



0.05 



Livetime 



Standard uncertainty arising from an 
estimated uncertainty of 0.007 % 
on the determination of the live time 



7- 10" 



Background 



Dilution factor 



Standard deviation on the determination of C^ 
determined via Monte Carlo simulation. A total of 
5 background data sets were constructed from random 
data arising from normal distributions defined by 
the average and standard deviation of experimental 
backgrounds at 4 efficiency points having 3 repetitions 
each; calculations were carried out with all 5 background 
data sets for a single experimental data set 

Uncertainty in C^ of solution GelAl due to uncertainty 
in gravimetrically-determined dilution factor between 
solutions in ampoules GelAl and GelAlD2 



5 • 10" 



0.04 



Combined (u^ = ^|llul) 



0.39 



Expanded {U^ = u^' k;k = 2) 



0.79 
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Table 6. Uncertainty components evaluated in the determination of the massic activity, C^, for ^^Ge solution^/ by liquid scintillation counting 
using CIEMAT/NIST H standard efficiency tracing method 



Component, w^- 



Comment 



Evaluation type 



Sample repeatability 



LS measurement 
reproducibility 

Mass determinations 



Dilution factor 

^^Ge decay corrections 
''Ge efficiency 
Livetime determinations 



Standard deviation of the mean on the 
determination of massic activity for a single LS 
cocktail (ji= iO determinations of C^ per source) 

Standard deviation on the determination of C^ for 
10 cocktails of 2 compositions 

Estimated standard uncertainty of ^^Ge mass for 
any single LS cocktail 

Uncertainty in C^ of solution GelAl due to 
uncertainty in gr avimetrically-det ermine d 
dilution factor between solutions in ampoules 
GelAl andGelAlD2 

Standard uncertainty in half-hfe (0.059 %) 
over the measurement decay interval 

Estimated uncertainty in C^ due to step size 
in CN2004 calculations 

Estimated uncertainty in the correction to the 
LS counting interval 



0.05 



0.24 



0.05 



0.04 



0.001 



0.65 



0.05 
(and PE)' 



Background 

Activity of ^H standard 

Branching ratios 



Estimated uncertainty due to an average 4 % 
uncertainty in background determination 

Estimated uncertainty due to 0.36 % uncertainty 
in ^H standard activity 

Estimated uncertainty due to uncertainty in 
branching ratios 



0.004 



0.18 



0.08 



Combined {Uq = ^Ym]) 



0.73 



Expanded (Uq = Uq' k; k = 2) 



1.45 



The relative uncertainty for this component is partially embodied (PE) in the relative standard uncertainties of the repeatability and reproducibility components. 
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Table 7. Uncertainty components evaluated in the determination of the massic activity, Q^, for Ge solution GelAl by ^-ray spectrometry using 
HPGe detectors 



Component, w^- 



Comment 



Evaluation type 



Measurement repeatability 

Efficiency curve 
Sample geometry 

Decay correction 
Decay data 

Dilution factor 



Standard deviation on determination of C^ 
for 3 repeated measurements of a single 
source at a single geometry 

Standard deviation of the mean on determination 
of detection efficiency for 4 sample geometries 

Typical uncertainty due to change of sample 
geometry (detector and source-to-detector distance) 
for a single counting source 

Standard uncertainty in half-Ufe (0.059 %) over 
the measurement decay interval 

Standard uncertainty (0.93 %) on emission 
probablity of 1078 keV gamma-ray in the 
decay of ^^Ge 

Estimated standard uncertainty in C^ of 
solution Gel Al due to uncertainty in 
gravimetrically-determined dilution factor 
between solutions in ampoules GelAl 
and GelAlD2 



B 
B 



0.98 

0.12 
0.33 

2.6 ■ 10-- 
0.93 

0.04 



Combined (Uq = ^1^u^) 



1.4 



Expanded (Uq = Uq' k; k = 2) 



2.8 



therefore account for all possible decay paths, 
including atomic rearrangements following electron 
capture. From a practical standpoint, a compromise 
between treating all possible paths and reasonable 
computation times must be made. While this certainly 
introduces some small amount of uncertainty, it is not 
expected that the weak contributions due to paths not 
considered in Fig. 4 would be significant, at least for 
the TDCR method. Instead, as seen in Table 5, the 
uncertainties on the input data play a very significant 
role. 

Because of the more direct nature of the measure- 
ment in the LTAC technique, the LTAC activity value 
for the solution in ^7 was adopted as the reference 
value for this study and was used in the calculation of 
the K-value for the NIST IC. The fact that the LTAC 
and TDCR measurements agree to within their respec- 
tive experimental uncertainties is encouraging, given 
the complexity of the TDCR efficiency calculation. 
Nonetheless, one would hope that improvements in the 
NIST TDCR spectrometer would lead to higher EC 



detection efficiencies, thereby providing better results 
in the measurement of radionuclides that decay by this 
mode. The CNET results indicate that some improve- 
ments in the method are still needed to be able to 
reliably measure nuclides that undergo EC decay. 

3.3 Determination of K- Value for NIST IC 

In order to avoid the need to perform a primary 
standardization every time a NIST-calibrated solution of 
^^Ge/^^Ga is required, we determined a calibration factor 
(K-value) for the NIST IC. This K-value is not to be 
confused with the coverage factor, k, applied to uncer- 
tainty evaluations. Using the LTAC reference activity 
value and the measured responses in the IC against 
radium reference sources (RR) 500B and 100, the 
K-values were found to be 2.695 x 10^ ± 1.7 x 10^ Bq 
and 5.032 x 10^ ± 3.2 x 10^ Bq, respectively The uncer- 
tainties on the K-values are expanded (k = 2) uncertainty 
and include relative standard uncertainty components 
due to the original primary standardization (0.29 %), 
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repeatability on 40 measurements in the IC (0.015 %), 
source mass (0.05 %), decay correction (0.002 %), and 
source positioning (0.1 %). 



4. Conclusion 

A solution containing ^^Ge in equilibrium with its 
decay daughter ^^Ga has been standardized for the first 
time at NIST, with a combined standard uncertainty of 
0.29 % using LTAC. Measurements made with two 
other LS techniques, TDCR and CNET, confirmed the 
LTAC result to within 1.1 % and 1.5 %, respectively. 
The differences between results obtained with the latter 
two methods and the LTAC technique indicate that 
improvements in the models and/or their applications 
are needed, particularly for EC nuclides. 

Data collected on the NIST 47C7 ionization chamber 
allowed for the determination of calibration factors for 
that chamber in the 5 mL NIST ampoule geometry, 
thereby enabling future calibrations of solutions having 
the same solution composition without the need for the 
measurements to be made by a primary method. 
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